PSO Method for Fitting Analytic Potential Energy Functions. Application to I-(H2O).

In this work a particle swarm optimization (PSO) algorithm was used to fit an analytic potential energy function to I-(H2O) intermolecular potential energy curves calculated with DFT/B97-1 theory. The analytic function is a sum of two-body terms, each written as a generalized sum of Buckingham and Lennard-Jones terms with only six parameters. Two models were used to describe the two-body terms between I- and H2O: a three-site model H2O and a four-site model including a ghost atom. The fits are compared with those obtained with a genetic/nonlinear least-squares algorithm. The ghost atom model significantly improves the fitting accuracy for both algorithms. The PSO fits are significantly more accurate and much less time-consuming than those obtained with the genetic/nonlinear least-squares algorithm. Eight I----H2O potential energy curves, fit with the PSO algorithm for the three- and four-site models, have RMSE of 1.37 and 0.22 kcal/mol and compute times of ∼20 and ∼68 min, respectively. The PSO fit for the four-site model is quite adequate for determining densities of states and partition functions for I-(H2O) n clusters at high energies and temperatures, respectively. The PSO algorithm was also applied to the eight potential energy curves, with the four-site model, for a short time ∼8 min fitting. The RMSE was small, only 0.37 kcal/mol, showing the high efficiency of the PSO algorithm with retention of a good fitting accuracy. The PSO algorithm is a good choice for fitting analytic potential energy functions, and for the work presented here was able to find an adequate fit to an I-(H2O) analytic intermolecular potential with a small number of parameters.

Electronic nature of zwitterionic alkali metal methanides, silanides and germanides - a combined experimental and computational approach.

Zwitterionic group 14 complexes of the alkali metals of formula [C(SiMe2OCH2CH2OMe)3M], (M-1), [Si(SiMe2OCH2CH2OMe)3M], (M-2), [Ge(SiMe2OCH2CH2OMe)3M], (M-3), where M = Li, Na or K, have been prepared, structurally characterized and their electronic nature was investigated by computational methods. Zwitterions M-2 and M-3 were synthesized via reactions of [Si(SiMe2OCH2CH2OMe)4] (2) and [Ge(SiMe2OCH2CH2OMe)4] (3) with MOBu t (M = Li, Na or K), resp., in almost quantitative yields, while M-1 were prepared from deprotonation of [HC(SiMe2OCH2CH2OMe)3] (1) with LiBu t , NaCH2Ph and KCH2Ph, resp. X-ray crystallographic studies and DFT calculations in the gas-phase, including calculations of the NPA charges confirm the zwitterionic nature of these compounds, with the alkali metal cations being rigidly locked and charge separated from the anion by the internal OCH2CH2OMe donor groups. Natural bond orbital (NBO) analysis and the second order perturbation theory analysis of the NBOs reveal significant hyperconjugative interactions in M-1-M-3, primarily between the lone pair and the antibonding Si-O orbitals, the extent of which decreases in the order M-1 > M-2 > M-3. The experimental basicities and the calculated gas-phase basicities of M-1-M-3 reveal the zwitterionic alkali metal methanides M-1 to be significantly stronger bases than the analogous silanides M-2 and germanium M-3.

Mechanistic details of energy transfer and soft landing in ala2-H(+) collisions with a F-SAM surface.

Previous chemical dynamics simulations (Phys. Chem. Chem. Phys., 2014, 16, 23769-23778) were analyzed to delineate atomistic details for collision of N-protonated dialanine (ala2-H(+)) with a C8 perfluorinated self-assembled monolayer (F-SAM) surface. Initial collision energies Ei of 5-70 eV and incident angles θi of 0° and 45°, with the surface normal, were considered. Four trajectory types were identified: (1) direct scattering; (2) temporary sticking/physisorption on top of the surface; (3) temporary penetration of the surface with additional physisorption on the surface; and (4) trapping on/in the surface, by physisorption or surface penetration, when the trajectory is terminated. Direct scattering increases from 12 to 100% as Ei is increased from 5 to 70 eV. For the direct scattering at 70 eV, at least one ala2-H(+) heavy atom penetrated the surface for all of the trajectories. For ∼33% of the trajectories all eleven of the ala2-H(+) heavy atoms penetrated the F-SAM at the time of deepest penetration. The importance of trapping decreased with increase in Ei, decreasing from 84 to 0% with Ei increase from 5 to 70 eV at θi = 0°. Somewhat surprisingly, the collisional energy transfers to the F-SAM surface and ala2-H(+) are overall insensitive to the trajectory type. The energy transfer to ala2-H(+) is primarily to vibration, with the transfer to rotation ∼10% or less. Adsorption and then trapping of ala2-H(+) is primarily a multi-step process, and the following five trapping mechanisms were identified: (i) physisorption-penetration-physisorption (phys-pen-phys); (ii) penetration-physisorption-penetration (pen-phys-pen); (iii) penetration-physisorption (pen-phys); (iv) physisorption-penetration (phys-pen); and (v) only physisorption (phys). For Ei = 5 eV, the pen-phys-pen, pen-phys, phys-pen, and phys trapping mechanisms have similar probabilities. For 13.5 eV, the phys-pen mechanism, important at 5 eV, is unimportant. The radius of gyration of ala2-H(+) was calculated once it is trapped on/in the F-SAM surface and trapping decreases the ion's compactness, in part by breaking hydrogen bonds. The ala2-H(+) + F-SAM simulations are compared with the penetration and trapping dynamics found in previous simulations of projectile + organic surface collisions.

Reaction dynamics of temperature-variable anion water clusters studied with crossed beams and by direct dynamics.

We present a study of the different product channels in the reactions of OH and OH-(H2O) with methyl iodide over a range of collision energies. Direct dynamics classical trajectory simulations are employed to obtain an atomistic comparison with the experimental results. For the experiments we have combined a crossed beam ion imaging setup with a multipole rf ion trap. The trap allows us to prepare the molecular and cluster ions with a controlled internal temperature and thus provides well-defined initial conditions for reaction experiments at low collision energy. Changing the internal temperature of the cluster ions was found to have a profound effect on their reactivity.

Higher-accuracy schemes for approximating the Hessian from electronic structure calculations in chemical dynamics simulations.

In this paper, we present a family of generally applicable schemes for updating the Hessian from electronic structure calculations based on an equation derived with compact finite difference (CFD). The CFD-based equation is of higher accuracy than the quasi-Newton equation on which existing generally applicable Hessian update schemes are based. Direct tests of Hessian update schemes, as well as dynamics simulations using an integrator incorporating Hessian update schemes, have shown four of the new schemes produce reliably higher accuracy than existing Hessian update schemes.

Imaging nucleophilic substitution dynamics.

Anion-molecule nucleophilic substitution (S(N)2) reactions are known for their rich reaction dynamics, caused by a complex potential energy surface with a submerged barrier and by weak coupling of the relevant rotational-vibrational quantum states. The dynamics of the S(N)2 reaction of Cl- + CH3I were uncovered in detail by using crossed molecular beam imaging. As a function of the collision energy, the transition from a complex-mediated reaction mechanism to direct backward scattering of the I- product was observed experimentally. Chemical dynamics calculations were performed that explain the observed energy transfer and reveal an indirect roundabout reaction mechanism involving CH3 rotation.

Trajectory studies of S(N)2 nucleophilic substitution. 8. Central barrier dynamics for gas phase Cl(-) + CH(3)Cl.

Quasiclassical direct dynamics trajectories, calculated at the MP2/6-31G level of theory, are used to study the central barrier dynamics for the C1(-) + CH(3)Cl S(N)2 reaction. Extensive recrossings of the central barrier are observed in the trajectories. The dynamics of the Cl(-)-CH(3)Cl complex is non-RRKM and transition state theory (TST) is predicted to be an inaccurate model for calculating the Cl(-) + CH(3)Cl S(N)2 rate constant. Direct dynamics trajectories also show that Cl(-) + CH(3)Cl trajectories, which collide backside along the S(N)2 reaction path, do not form the Cl(-)-CH(3)Cl complex. This arises from weak coupling between the Cl(-)-CH(3)Cl intermolecular and CH(3)Cl intramolecular modes. The trajectory results are very similar to those of a previous trajectory study, based on a HF/6-31G* analytic potential energy function, which gives a less accurate representation of the central barrier region of the Cl(-) + CH(3)Cl reaction than does the MP2/6-31G* level of theory used here. Experiments are suggested for investigating the non-RRKM and non-TST dynamics predicted by the trajectories.

Simulations of Gas-Phase Chemical Reactions: Applications to SN2 Nucleophilic Substitution.

Computer simulations and animations of the motion of atoms as a chemical reaction proceeds give a detailed picture of how the reaction occurs at a microscopic level. This information is particularly useful for testing the accuracy of statistical models, which are used to calculate various attributes of chemical reactions. Such simulations and animations, in concert with experimental and ab initio studies, have begun to provide a microscopic picture of the intimate details of a particular class of gas-phase ion-molecule bimolecular reactions known as S(N)2 nucleophilic substitution. In these reactions, a nucleophile is displaced from a molecule by another nucleophile. The dynamical model of S(N)2 reactions that emerges from the computer studies, and its relation to statistical theories, is discussed here.